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ABSTRACT 


A  general  method  to  compute  the  parameters  of  microstrip 
transmission  lines  is  described.  Here,  the  values  of  the 
potential  function  are  estimated  by  solving  the  Laplace 
equation  using  the  Monte  Carlo  method.  The  electric  field 
strength  is  then  computed  by  differentiating  the  potential 
function.  From  the  electric  field  distribution,  the  charge 
density  on  the  surface  of  center  conductor,  the  associated 
capacitance  per  unit  length  of  the  line,  the  characteristic 
impedance,  and  the  conductor  loss  of  the  line  is  computed. 

A  Fortran  program  has  been  prepared  and  the  flow  chart  of 
the  program  is  included  and  described  in  this  paper.  The 
program  was  used  to  estimate  the  characteristic  impedance 
of  lines  of  three  different  geometries.  Results  are  within 
ten  percent  of  the  experimentally  measured  values. 


iii 
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INTRODUCTION 

This  paper  is  concerned  with  a  general  method  to  compute  the 
characteristic  impedance  and  conductor  loss  of  microstrip  transmission 
lines  consisting  of  either  a  conductor  atrip  above  an  infinite  ground 
plane  or  a  strip  between  two  ground  planes.  It  has  been  proven,  both 
experimentally  and  theoretically,  that  for  all  practical  purpose  we  may 
assume  that  the  conductors  are  perfect  and  are  imbedded  in  homogeneous 
and  lossless  dielectric  material*  Therefore,  TEM  mode  may  propagate  in 
these  types  of  lines.  Thus  the  analysis  of  the  parameters  of  these 
lines  is  essentially  reduced  to  the  investigation  of  the  electric  field 
distribution  in  the  transverse  cross-sectional  plane  of  the  line.  From 
the  field  distribution,  the  charge  density  on  the  surface  of  the  con- 
d-ictoi.  capacitance  per  unit  length  of  the  line,  and  the  characte*  • 
istic  impe.  ancr  tr*  line  can  be  determined.  This  problem  is,  there¬ 
fore  one  :>f  static  elec;  1  ^magnetism. 

However,  even  after  this  simplification  is  n-:.da,  to  derive  an 
analytical  expi  ^ou  io.  the  parameters  of  such  micrortrip  liner  is 
still  rather  difficult.  By  making  approximations  in  one  way  or  another, 
many  people  have  worked  out  solutions  which  give  the  parameters  of  lines 
with  specific  geomerrfcs.  For  the  lines  where  the  widt  of  the  conductor 
strip  (hereon  denoted  as  ’’*)  is  large  comparing  with  the  spacing  hecween 
the  strip  and  the  ground  plane  (hereon  denoted  as  H),  e.g,,  B  Jfc,  i.O  9 , 
the  edge  effect  is  negligible  at  the  center  of  the  strip.  This  case  has 
been  solved  by  using  the  conformal  mapping  method.  However,  in  the  cases 
where  B  is  not  large  compared  with  H,  the  solution  thus  obtained  are 
certainly  not  accurate.  Ihij  problem  has  also  been  solved  using  conformal 
mapping  method  for  any  ratio  of  B  and  H;  however  in  these  instances,  the 
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thickness  of  the  conductor  D  is  assumed  to  be  negligibly  smell  compered 
with  other  dimensions. 

The  anelytlcal  soltlon  for  the  most  general  case  is  just  too 
complicated,  especially  since  we  want  this  solution  to  be  valid  for 
both  cases  of  open  strip  lines  and  3&ndwlch  lines.  In  view  of  the  other 
approximation  methods,  it  seems  to  be  most  suitable  to  solve  our  problem 
by  solving  the  two  dimensional  Laplace’s  equation  numerically  since  a 
computer  program  may  be  prepared  for  the  most  general  case  so  as  to 
include  all  possible  geometries  of  the  line. 

Theoretically,  in  our  problem,  the  boundaries  includes  points 

at  infinity.  Because  only  at  these  p  ints  do  the  flalda  due  to  the 

charges  on  tne  conductor  strip  vanish.  However,  at  a  distance  large 

enough  away  from  the  strip,  the  electr'c  tieJU  is  practically  zero. 

Therefore,  we  can  approximate  the  boun. *ry  conditions  as  that  shown  in 

Fig.  1,  where  and  p  are  large  numbrrs  so  c'  osen  that  at  y  »  oi  and 

x  «  +  f  the  field  strength  is  practically  /.ero.  Also,  since  all 

2H 

parameters  we  want  to  compute  except  t:»*  pcwev  loss  per  unit  length  of 
the  line,  are  functions  of  the  ratios  B/H  and  bf H,  no  generality  is 
lost  by  letting  H  *  1.  Then,  f:r  any  given  the  power  loss  per  unit 

length  can  be  easily  obtained  b>  multiplying  the  power  loss  obtained 
from  this  geometry  by  the  proper  multiplication  factor. 

i  mrfaoe  charge  de.ui  v  t  a  center  conductor  and  the 
associated  capacitance  can  be  determined  from  the  field  distribution 
just  outside  of  the  strip.  From  the  value  of  the  capacitance  and  the 
charge  density,  the  characteristic  impedance  and  power  flow  may  be 
determined.  Therefore,  it  is  only  necessary  to  find  the  electric 
potential  at  points  cios^:  tc  the  cen*.  tx  conductor.  For  this  reason,  the 
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Monts  Csrlo  method  was  chosen  to  solve  the  Laplace  equation,  because 
this  method  enables  us  to  obtain  the  solution  of  the  Laplace's  equatlc 
at  these  points  without  having  to  obtain  the  solution  at  any  other  points 
at  the  same  time. 


Fig.  I 
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Fig.  2 


DESCRIPTION  OF  MONTE  CARLO  METHOD  FOR  SOLVING  THE  IAPIACE’S  EQUATION 

It  may  be  helpful  to  describe  the  well  known  Monte  Carlo 
Method  briefly.  Given  a  Poisson5 s  equation 


+ 


_^v 

iz2 


F  (x,  y,  z) 


where  F  (x,  y,  z)  is  e  function  specified  within  e  closed  boundary 
surface  S,  V(x,  y,  z)  is  the  potential  at  the  point  whose  coordinates 
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are  x,  y ,  z.  On  the  boundary  surface  S,  the  value  of  the  potential 
function,  V,  is  given. 

Consider  that  the  region  inside  the  surface,  S,  is  divided 
up  into  small  cubical  lattices  of  sides  h.  And  the  boundary  surface 
is  replaced  by  a  set  of  boundary  points  J-  as  shown  in  Fig.  2. 

This  set  of  points,  |  includes  all>nearest  exterior  neighbors 

to  all  lattice  points  witnin  S  and  the  nearest  exterior  points 
diagonally  across  the  lattice  points  in  S.  The  valuss  of  the  potential 
function  V  at  all  boundary  points  j^P^  j  are  given  as-^^,  y^,  z^) 

Then  the  poisson's  equation  may  be  approximated  as  t.ne  difference  equatior 


A  V  +  AV+AV-F  (x.y.z)  . 
x  y  z 


(1) 


where 


AXV  "  ^2  {V  *x  +  +  v  (x-h,y,z*  -  2tf  <x,y. 

Ayv  "  {v  (X,  y  +  h,  z>  +  V  (X,  Vh,  z.)  -  2V  (x,  y,  z4 

{V  (x*  v»  *  +  h)  +  V  (x,  y,  z*h)  -  2V  (x,  y,  z.>j 


Denote  the  point  whose  coordinates  are  x,  y,  z  as  P  and  its  closest 
neighboring  points  as  P^  ?2  . Pfi.  Equation  (1)  becomes: 


V  (p) 


(2) 


Note  that  the  difference  equation  (2)  gives  the  relation  between  the 
potential  at  any  point  P  and  those  at  its  adjacent  points.  In  the 
Monte  Carlo  Method,  we  shall  consider  the  coefficients  of  the  difference 


equation  as  transition  probabilities  from  one  point  P  to  its  six 
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closest  neighbors.  Then  for  a  set  of  fictitious  particles  taking 
random  steps  at  the  point  P,  with  the  directions  of  their  steps 
governed  by  random  sampling  technique,  the  probability  for  these 
particles  to  step  to  any  one  of  the  neighboring  points  is  given  by 
the  coefficient  of  the  difference  equation.  In  our  case,  the  par¬ 
ticles  have  equal  probability  to  step  to  any  one  of  the  six  nearest 
neighboring  points. 

Now  consider  a  particle  is  allowed  to  take  X  random  walks 

where  all  walks  starts  at  a  lattice  point  Pq  and  terminates  whenever 

the  particle  reaches  the  boundary.  As  the  particle  is  momentarily 

at  a  point  P,  the  probability  that  it  will  step  to  each  one  of  the 

six  closest  neighbors  of  P  is  1/6.  A  tally  of  the  value  of  the 

potential  function  made  for  each  walk  will  depend  on  the  transition 

probability  between  successive  points  and  the  value  of  F(P)  at  the 

points  which  the  particle  reached  in  the  walk.  It  also  depends  on 

the  value  of  the  potential  at  the  boundary  point, y^,  z^), where 

th 

the  walk  Is  terminated.  Here  for  the  j  walk,  the  expected  value  of 
the  tally  Z^  is  given  as 

Zj  -  -  h2-  F  (P±)  +  4>(Xj.  Yj »  Zj)  .  (3) 

where  (Xj ,  y  ,  z^)  is  the  boundary  point  at  which  the  j**1  walk  is 
terminated  and  P^’s  are  points  reached  in  the  walk.  The  solution  of 
equation  (2)  at  the  point  P.  can  be  estimated  as: 

5  <"„>  -  i  I  '] 

j-l 


\ 
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And  the  error  of  the  solution  1.  estimated  by  the  variance 
c2  (V)  •  £  (2j  -  V2) 

It  it  obvious  here  that  the  larger  X  is.  the  higher  the  accuracy  will 
be. 

As  mentioned  in  the  introduction,  the  analysis  of  the  micro¬ 
strip  line  is  simplified  as  a  two-dimensional  problem.  Also,  with 
the  dielectric  being  homogeneous  and  lossless  and  the  conductors  being 
perfect,  the  function  F(p)  in  the  poisson's  aquation  is  sero  at  all 
points  in  the  boundary.  Hence  aquations  (2)  and  (3)  become: 


V(P)'J  ^  V  (Pj)  .  . .  (4) 

And 

Zj  "  ^  <XJ*  yj>  .  (5) 


DESCRIPTION  OF  THE  COMPUTATION  PROCEDURE  AND  FLOW  CHART 

First  let  us  redraw  the  cross-section  ot  a  microstrip  line 

in  Fig.  3,  where  all  dimensions  all  indicated  bv  the  numher  of  meshes 

H 

of  aidaa  h.  It  la  clear  tnat  h  is  equal  to  —  .  The:  coordinates  of 

N 

each  of  the  boundary  points  are  Indicated  in  this  figure  in  terms  of 
the  indices  I  and  J  which  denote  the  lattice  point  representing  the 
boundary  according  to  the  rule  described.  Obviously,  from  the  rule, 

(  )  is  the  largest  integer  &  |  N.  ^11^1  and  Ara 

largest  Integers  aqual  to  or  smaller  than  jjg  and  -g-  respectively. 

BN 

With  ^  being  large'  and  arbitrary,  we  can  always  make  |J  —  an  integer. 


MfW+! _ _ 

folk  SHMPWICH  UMfc&i 
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It- nay  «1m  Ix  helpful  to  list  here  •  few  well  known-  rel-atlone 
which  give  us  the  parameters  of  interest  in  terms  of  the  scaler  potential, 
V. 

The  electric  field  intensity,  E  ■  -  Grad  V  . 

The  charge  density  . .  D  ■  4  E  . 

The  total  surface  charge  on  the  strip  per  unit  length,  q , 

■JaCDFA  (D*n)ds  .  (c) 

(where  t  is  a  unit  vector  pointing  outward  of  and 
being  perpendicular  to  the  surface  of  center  conductor .) 


Therefore,  q  oc 
The  characteristic  impedance, 


<d) 


Alpha  *  max.  possible  power  loss/power  transmitted 


From  these  relations  we  can  derive  all  parameters  of  interest  once 
the  electric  potential  at  points  just  outside  of  the  center  strip  is 
known.  Because  of  the  symmetry  of  the  geometry  of  the  line,  it  is 
sufficient  to  compute  the  potential  at  points  whose  coordinates  are: 
1)  from  I  -  LA  to  I  -  LAP  and  J  -  N,  2)  from  I  -  LA  to  I  -  LAP  and 
J  -  NDI  and  3)i-  LAP+1  and  J  «  N  to  J  -  ND1. 

We  shall  describe  the  flow  chart  in  datail. 
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Box  1. 


Box  2. 


Box  3 


Box  4 


Box  5 


Box  6 


Input  data: 

B  D  V 

a)  Dimensional  constants:  H  (in  meters),  — ,  -,  N,  and 

(for  sandwich  lines  only) 

b)  Number  of  walks  to  be  made,  X 

c)  Arbitrarily  chosen  constants  o(  (for  open  line)  and  ^ 

d)  Dielectric  constants,  &  ,  permeability  fx ,  conductivity 
of  the  conductors  ce,  and  frequency  at  which  the  power 
loss  is  of  interest, 

a)  Compute  constants  and  boundary  points  LA,  NP,  L,  etc. 
and  express  them  in  fixed  point  numbers  conforming  with 
the  rule  mentioned  above. 

Register  Y  is  set  to  one,  and  solution  point  P  (I,J  )  is 

assigned  to  be  P  (LA,N) .  Register,  P,  for  integral  of 

electric  field  strength  E  and  Register, Q,  for  integral  of 

2 

the  square  of  the  electric  field  strength,  E  ,  are  set  to 
be  zero. 

Sum  register,  G,  for  tally  and  sum  of  square  register  C 
2 

for  tally  are  set  to  zero.  The  counterIX  indexing  the 
walks  is  set  to  one. 

Current  point  counters  and  Jg  are  set  to  the  indices  of 
the  solution  point  P(I,J  ).  Register  Z  for  tally  during  a 
walk  is  set  to  zero. 

Two  random  number  U#  and  are  generated.  Here  the  sub¬ 
routine  RANF  already  in  Fortran  Library  which  generates 
floating  random  numbers  of  magnitude  between  0  and  1  with 


even  distribution  are  used. 
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Box  7.  impute  h  and  A  J#  according  to  the  table  below: 


>>e 


£ 

s 


I 


s 


J 


0.5 

>  0.5 

1 

0 

0.  5 

<.  0.5 

0 

1 

0.5 

>  0,5 

0 

-1 

0.5 

<  0.5 

-1 

0 

Box  8. 


The  current  point  are  advanced  to  where  I  »  Ic  +  A  1^, 


J  +  A  J 

s  s 


liox  9. 


Box  10. 


Box  11. 


Box  12. 


Ecx  13. 


Box  15. 


Test,  if  the  new  current  point  coincides  with  any  of  the 
bo  ndery  points.  Tf  it  does  not,  step  back  to  Box  £.  If 
it  does,  proceed  to  Box  10. 

The  boundary  value  at  the  boundary  point  reached  is  trans¬ 
ferred  to  the  tally  register  Z. 

The  content  cf  the  tally  register  Z  is  added  tc  that  of  ere 
sum  register  G.  And  the  square  of  the  content  of  register 
Z  is  added  to  that  of  the  sum  of  square  reg:» C. 

If  ll  walks  ire  completed,  proceed  co  Bon.  13.  Otherwise, 
srep  to  Box  14  where  the  register IX  is  increased  by  one. 
Thu...  'up  baci«  t  *  S^x  5  co  ct*rt  a  walk, 

Esti»r!*e  of  the  solution  at  the  point  (I3  J)  is  computed  ard 
store  in  V(I ,  J) . 

The  'criance  of  the  estimated  solution  is  computed. 
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Box  16.  If  the  content  of  the  register  Y  is  one  (solution  points 
are  on  line  E'D*  in  Fig.  3,  proceed  to  Box  17.  If  it  is 
two  (solution  points  are  on  line  B‘C')  go  to  Box  26.  If 
it  is  3  (solution  in  points  are  on  line  C'D*  )  go  to  box  27. 

Box  17.  Compute  the  electric  field  strength  using  the  approximated 
equation  E  (I,J  )  ■  Z  v(*i  OUl  #  Here  (on  the  line 

E'D')  the  x-directional  electric  field  is  zero.  V(I,  J+l) 
is  given  as  boundary  value. 

Box  18.  If  I  is  equal  to  LA,  increase  the  register  I  by  one  and  step 
back  to  Box  4  to  compute  the  solution  for  the  next  point. 
Otherwise,  with  solutions  for  more  than  one  point  are  avail- 
able,  go  to  Box  19. 

Box  19.  Compute  the  average  field  strength  between  the  two  adjacent 
points  (I,  J)  and  (1-1,  J). 


Box  20.  Compute  the  line  integral 


Box  21. 


if  |E| 

J acdfa 


ds  by  using  the 


' ACDFA 

trapezoidal  rule  of  integration. 


4 

J  AC 


Compute  the  line  integral  I  |E(  ds  by  using  the 

/acdfa 


trapezoidal  rule  of  integration. 

Box  22.  If  the  solutions  for  all  points  with  coordinates  I  -  LA  to 

I  ■  LAP  and  J  ■  N  are  obtained,  proceed  to  Box  23.  If  other¬ 
wise,  stop  back  to  Box  18b. 

Box  23.  If  register  Y  contains  the  number  one,  go  to  Box  24.  If  it 
contains  the  number  two,  go  to  Box  25. 

Box  24.  Solution  point  registers  I  and  J  are  assigned  to  be  LA  and 
ND1.  And  the  register  Y  is  increased  by  one.  Then,  step 


back  to  Box  4. 
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•  2ti  ■  v*m«»ite  electric  field  strength  at  points  on  the  line  B'C' 
in  I..,  .  3.  Here  the  point  (I,  J-l)  is  on  the  surface  of 
"he  center  strip. 

Box  27.  On:..- at  a  electric  field  strength  at  point  on  line  C  'D'  where 

the  vertical  electric  field  is  zero  (1*1,  J>  are  boundary 

points . 

Boa  2.-  If  J  is  smaller  than  NP+1  (the  solution  for  only  one  point 
on  the  line  C'D*  is  available),  step  back  to  Box  28b  where 
the  register  J  is  increased  by  one.  Then,  step  back  to  Box  4 
to  compute  the  solution  of  another  point.  If  otherwise, 
proceed  to  Box  29. 

Boxes'29,  30  &  31.  Same  as  boxas  19,  20,  and  21  respectively. 

Box  32.  It  the  solution  of  all  points  are  obtained  proceed  to  Box  33. 
If  net,  step  back  to  28b. 

Box  33.  Compute  the  characteristic  impedance  of  the  line. 

3*.  Co:r.puc«  A. 

Bex  35.  Vu.nt  out  the  estimated  potential  at  all  points  on  the  lines 
B'C'.  C'D'  and  D'E-  with  the  variances  of  the  estimated 
"it  ions.  The  characteristic  Impedance  of  the  line  and  the 
:  Uo  .LPHA. 


LIMITATION  vc  THE  PROGRAM 

Du.  t o  the  limitation  of  the  storage  space,  using  an 
IBM  7090  cciputer  with  32K  words  storage,  the  maximum  number  of  meshes 


W-5490 


1 6 


the  dimensions  H,  B/H  end  D/H  may  be  divided  into  ere  ee  follow*: 


otN  4  2,000 

f  4  6,000 

Therefore  H  ■  2>000/oL  ,  B/H  «  6,000/(jbN) 

Thle  program  was  used  to  estimate  the  characteristic 
impedance  of  open  strip  lines  of  the  dimansions:  B/H  *  1,0,  2,0, 
and  3,0  and  D/H  ■  0,3*  H  was  divided  into  ten  increments.  Both  cl 
and  ji  were  chosen  to  be  ten.  For  one  hundred  random  walks,  the 
results  obtained  are  within  10%  of  the  experimentally  measured  values. 
Although  dividing  the  region  into  fewer  meshes  and  increasing  the 
number  of  walks  will  certainly  Improve  the  accuracy,  the  computing 
time  have  been  found  to  be  extremely  high  (say,  one  hour  on  the  IBM 
7090). 

More  accurate  results  will  be  obtained  in  a  reasonable  amount 
of  computing  time  in  the  cases  where  the  lines  are  the  sandwich  type. 
The  reason  is  that  with  the  spacing  between  the  two  horizontal  boundary 
lines  at  ground  potential  being  not  more  than  two  to  five  times  the 
spacing  H;  the  mean  distance  of  the  random  walks  is  shorter;  therefore 
each  walk  will  take  less  time. 


<0. 


Jane  tf.  S.  Liu 
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